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ABSTRACT 

We constructed some main-sequence mergers from case A binary evolution and stud- 
ied their characteristics via Eggleton's stellar evolution code. Both total mass and 
orbital angular momentum are conservative in our binary evolutions. Assuming that 
the matter from the secondary homogeneously mixes with the envelope of the primary 
and that no mass are lost from the system during the merger process, we found that 
some mergers might be on the left of the zero-age main sequence as defined by normal 
surface composition (i.e helium content Y = 0.28 with metallicity Z — 0.02 for Pop I) 
on a colour-magnitude diagram(CMD) because of enhanced surface helium content. 
The study also shows that central hydrogen content of the mergers is independent of 
mass. Our simple models provide a possible way to explain a few blue stragglers (BSs) 
observed on the left of zero-age main sequence in some clusters, but the concentration 
toward the blue side of the main sequence with decreasing mass predicted by Sandquist 
et al. will not appear in our models. The products with little central hydrogen in our 
models are probably subgiants when they are formed, since the primaries in the pro- 
genitors also have little central hydrogen and will likely leave the main sequence during 
merger process. As a consequence, we fit the formula of magnitude and B — V of 
the mergers when they return back to thermal equilibrium with maximum error 0.29 
and 0.037, respectively. 

Employing the consequences above, we performed Monte Carlo simulations to 
examine our models in an old open cluster NGC 2682 and an intermediate-age cluster 
NGC 2660. Angular momentum loss (AML) of low mass binaries is very important in 
NGC 2682 and its effect was estimated in a simple way. In NGC 2682, binary mergers 
from our models cover the region with high luminosity and those from AML are located 
in the region with low luminosity, existing a certain width. The BSs from AML are 
much more than those from our models, indicating that AML of low mass binaries 
makes a major contribution to BSs in this old cluster. Our models are corresponding 
for several BSs in NGC 2660. At the region with the most opportunity on the CMD, 
however, no BSs have been observed at present. Our results are well-matched to 
the observations if there is 0.5Mq of mass loss in the merger process, but 
a physical mechanism for this much mass loss is a problem. 
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1 INTRODUCTION 

Much evidence shows that primordial binaries make 
an important contribution to blue stragglers (BSs) 
(IFerrar o et al. 20031 [Davies, Piotto fc Angeli 2004| 

[Mapclli et al.2004p . At present, a few BSs, i.e. F190, 
6 Car, have already been confirmed to be in binaries by 
observations, and their formation may be interpreted by 
mass transfer between the components of a binary. Whereas 
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in intermediate-age and old open and globular clusters, the 
number of observed close binaries among well-studied BSs 
is consistent with the hypothesis of binary coalescence. For 
example, Mateo et al. (|1990| made a comparison of the 
number of close binaries with the total number of BSs in 
NGC 5466 and found that it is an acceptable claim that all 
non-eclipsing BSs are formed as the result of mergers of the 
components in close binaries, though the possibility of other 
mechanisms to produce BSs cannot be ruled out due to the 
large uncertainties in their analysis. Monte-Carlo simula- 
tions of binary stellar evolution l|Pols fc Marinus 1994P also 



2 Chen and Han 



show that binary coalescence may be an important channel 
to form BSs in some clusters (e.g. with an age greater than 
40 Myr). Meanwhile, the arguments in theory show that W 
UMa binaries (low-mass contact binaries) must eventually 
merge into a single star HWebbink 19761 IWebbink 19851 
ITutukov fc Yungelson 19871 IMateo et al.l990p . Observa- 
tionally, the lack of radial velocity variations for most 
BSs further indicates that binary coalescence may be 
more important than mass transfer for BS formation 
UStryker 1993] IPols fc Marinus 1994|) . FK Comae stars 
are generally considered to be direct evidence for binary 
coalescence ( [Stryker 1993p . The smallest mass ratio of 
components among observed W UMa systems to date is 
about 0.06. All of the above show that it is important to 
study the remnants of close binaries. However the merge 
process is complicated and the physics during the process 
is still uncertain. Recently, Andronov, Pinsonneault fc 
Terndrup l|2006p studied the mergers of close primordial 
binaries by employing the angular momentum loss rate 
inferred from the spindown of open cluster stars. Their 
study shows that main sequence mergers can account for 
the observed number of single BSs in M67 and that such 
mergers are responsible for at least one third of the BSs 
in open clusters older than 1 Gyr. The physics of mergers 
are limiting case treatments in the study of Andronov, 
Pinsonneault fc Terndrup (|2006|l . Based on previous studies 
of contact binaries and some assumptions, we construct a 
series of merger models in this paper, to study the structure 
and evolution of the models and show some comparisons 
with observations. 



Case A binary evolution has been well studied by Nel- 
son & Eggleton (|200ip . They defined six major subtypes 
for the evolution (AD, AR, AS, AE, AL and AN) and 
two rare cases (AG and AB). Three of the subtypes (AD, 
AR, AS) lead the binary contact as both components are 
main-sequence stars and two cases (AE and AG) reach con- 
tact with one or both components having left the main se- 
quence. As there is no description for weird objects except 
for two merged main-sequence stars, merger products (ex- 
cept for two main-sequences stars) are generally assumed 
to have terminated their evolution l|Pols fc Marinus 1994p . 
i.e. they have left the main sequence and cannot be rec- 
ognized as BSs. Here we are interested in the cases of two 
main-sequence stars, i.e. cases AS, AR and AD. If tdyn, ^kh, 
tMS represent the dynamic timescale, thermal timescale and 
main sequence timescale of the primary (the initial massive 
star,*l), respectively, the following shows a simple defini- 
tion of the three evolutionary cases: AD-dynamic Roche 
lobe overflow (RLOF), M > M/Uyn; AR-rapid evolution 
to contact, M > M/tKH,t cont act -tRLOF(*l) < O.Ums(*1); 
AS-slow evolution to contact, tcontact — iRLOF(*l) > 
0.1tMs(*l), where tRLOF and tcontact are the ages at which 
RLOF begins and the binary comes into contact, respec- 
tively. In case AD, the core of the secondary spirals in 
quickly and stays in the center of the merger. The merger 
then has a chemical composition similar to that of the 
primary, resembling the result of smoothed particle hy- 



drodynamic calculations ( 


Lombardi, Rasio & Shapirol996| 
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s et al. 2001 p. We therefore stud- 



ied just the systems in cases AR and AS for this work. 



2 ASSUMPTIONS 

Using the stellar evolution code devised by Eggleton 
(fWn which has been updated with the lat- 

est physics over the last three decades (IHan et al. 19941 
IPols et al. 19951 IPols et al. 1998) . we re-calculate the mod- 
els of cases AS and AR with primary masses between 0.89 
and 2Mq until the systems become contact binaries. The 
structures of the primaries and the compositions of the sec- 
ondaries are stored to construct the merger remnants. 

Before the system comes into contact, the accreting 
matter is assumed to be deposited onto the surface of the 
secondary with zero falling velocity and distributed homo- 
geneously all over the outer layers. The change of chemical 
composition on the secondary's surface caused by the ac- 
creting matter is 

dX,/dt = {dM/dt)/[{dM/dt)dt + Ms] ■ (X,a - X,s), (1) 

where dM/dt is the mass accretion rate, Xia and Xi^ are 
element abundances of the accreting matter and of the sec- 
ondary's surface for species i, respectively, and Ms is the 
mass of the outermost layer of the secondary. The value 
of Ms will change with the moving of the non-Lagrangian 
mesh as well as the chosen model resolution, but it is so 
smaU (~ 10"" - 10"^^Mq) in comparison with [dM / dt)At 
(~ 10"^ - IO-'^Mq) during RLOF that we may ignore the 
effect of various Ms on element abundances. Before and af- 
ter RLOF, we get dXi/dt — from the equation, which is 
reasonable in the absence of mixing HChen fc Han 2004|l . 

The merger models are constructed based on the fol- 
lowing assumptions: (i) contact binaries with two main- 
sequence components coalesce finally and the changes of 
structures of individual components during coalescence are 
ignored; (ii) the matter of the secondary is homogeneously 
mixed with that of the primary beyond the core-envelope 
transition point, which separates the core and the envelope 
of the mass donor; (iii) the system mass is conserved. 

Firstly, we present a brief discussion on these assump- 
tions. Webbink H1976P studied the evolutionary fate of low- 
mass contact binaries, and found that a system cannot sus- 
tain its binary character beyond the limits set by marginal 
contact evolution (/x = Mi/{Mi+M2) = 1.0). He stated that 
a contact binary will very likely coalesce as the primary is 
still on the main sequence in a real system. Up to now, it is 
widely believed that case AD probably leads to common en- 
velope, spiral-in, and coalescence on quite a short timescale. 
The final consequences of AS and AR are not very clear, but 
Eggleton ( 2000 ) pointed out that systems undergoing AR or 
AS evolution may maintain a shallow contact (perhaps in- 
termittently) as the mass ratio becomes more extreme, and 
finally coalesce. Recent study on W UMa (Li, Zhang & Han, 
2005) also shows that these systems will be eventually coa- 
lescence. The merged timescale, i.e. the time from a binary 
contact to coalescence, is important here. If it is too long, 
the structures of both components will change remarkably 
and the system may have not completed coalescence within 
the cluster age. There are many conflicting estimates for 
the timescale, however, from observations and theoretical 
models of the merger process. Early observational es- 
timates range from 10^-10* yr in various environments 
Hvan't Veer 19791 [Eggen fc Iben 1989] ). Th e following study 
explored the average age about 5 x 10* yr (|van't Veer 19941 
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IDryomova fc Svechnikov 2002) . Bilir et al. (|2005|l pointed 
out that the age difference between fieid contact binaries 
and chromosphericaiiy active binaries, 1.61 Gyr, is likely an 
upper limit for the contact stage by assuming an equilib- 
rium in the Galaxy, whereas the study of W UMa by Li, 
Han & Zhang 1) 2004p suggested a much longer timescale, 
about 7 Gyr. We adopt the empirically estimated values in 
this paper (i.e from 5 x 10^ to 1 x 10® yrs) and ignore 
the changes of structure of individual components during 
merger process. For low-mass contact binaries, the common 
envelope is convective HWebbink 1976^ . and the matter in 
it is thus homogeneous. If a system mimics shallow contact 
during coalescence, it is reasonable to assume that the mat- 
ter of the secondary mixes with the envelope homogeneously. 
Van't Veer (|1997p found that the mass loss from the system 
during coalescence is at a rate of about 2 x lO~^°M0yr~^ 
by observations. If we consider that the coalescence time 
is 5 X 10* yr in a binary, only O.IM0 is lost from the sys- 
tem as the binary finally becomes a single star. We then 
roughly assume that the mass is conservative during coa- 
lescence. However mass loss might be an important way to 
carry orbital angular momentum away from the binary in 
this process. 

Secondly, we discuss the choice of the core-envelope 
transition point which separates the core and the envelope 
in the primary. Many characteristics of the merger are rele- 
vant to the choice, e.g. the chemical composition in the en- 
velope, evolutionary track on Hertzsprung-Russel diagram, 
and some observational characteristics. Unfortunately, one 
cannot find the core-envelope transition point in a main- 
sequence star as easily as in evolved stars because the den- 
sity profile, as well as many other thermodynamic quantities 
(entropy, pressure, temperature etc.), is smooth and does 
not have a deep gradient for main-sequence stars. Chen & 
Han (|2005p studied the influences of core-envelope transition 
point on the mergers of contact binaries with two main- 
sequence components. They found that one may ignore the 
effects which result from different choices of the transition 
point on colours and magnitudes of the merger if it is outside 
the nuclear reaction region of the primary, which is com- 
monly considered as the nearest boundary of the secondary 
reaches in cases AS and AR. In this paper, the core-envelope 
transition is determined as the point within which the core 
produces 99 per cent of total luminosity. This choice is gener- 
ally outside the nuclear reaction regions and has little effect 
on the flnal results. 

Finally the merger remnant is constructed as follows: it 
has the total mass of the system and a chemical composition 
within Mic similar to the core of the primary. The chemical 
composition in the envelope of the merger is given by 



X, = (M,2 4- Mab)/(M2 + mt). 



(2) 



where Mi2 and Mub denote the total masses of species 
i of the secondary and of the primary's envelope, respec- 
tively. TTib is the envelope mass of the primary. There 
might be a region in which the helium abundance is less 
than that of the outer region. The matter in this region 
then has a lower mean molecular weight than that in the 
outer region. This results in secular instability and thermo- 
haline mixing ( [Kippenhanhn, Ruschenplatt fc Thomasl980| 
lUlrich 1972|) . We include it as a diffusion process in our code 
(|Chen fc Han 2004|) . 



In the models of Nelson & Eggleton l|200ip . both to- 
tal mass and angular momentum are conservative. It was 
mentioned by the authors, however, that these assumptions 
were only reasonable for a restricted range of intermediate 
masses, i.e spectra from about GO to Bl and luminosity 
class III-V. Observationally, some low mass binaries with 
late-type components show clear signs of magnetic activity, 
which indicates that the systems evolve by way of a scenario 
implying angular momentum loss (AML) by magnetic brak- 
ing l|Mestel 1984[l . Magnetized stellar winds probably do not 
carry off much mass, but they are rich in angular momen- 
tum because of magnetic linkage to the binaries. For close 
binaries, rotation is expected to synchronize with orbital 
period, so AML is at the expense of the orbital angular mo- 
mentum, resulting in orbital decaying and the components 
approaching each other. A detached binaries, then, may be- 
come contact and flnally coalesce at or before the cluster age 
( Stc pien 1995[ ). There are a number of subjects including the 
treatment of AML (|Li, Han fc Zhang 2004| [Stepien 200"6| 



IMicheal fc Kevin 20061 IDemircan et al. 2006|) . For simplic- 
ity, the conservative assumption is also adopted in our bi- 
nary evolutions. In old clusters, however, AML of low mass 
binaries is very important and we estimate its impor- 
tance in another way (see section 4.2). 



3 EVOLUTIONARY RESULTS 

A set of binaries undergoing AS and AR evolution from 
Nelson fc Eggleton (|200ip are choosen to study the charac- 
teristics of the merger products and their connections with 
blue stragglers. Table Al gives the initial parameters of 
the binary systems, their RLOF information and the struc- 
tures and evolutionary consequences of the mergers. The 
first three columns contain the initial mass of the primary 
Mii, the initial mass ratio gi (the primary to the secondary) 
and the initial orbital period Pi in logarithmic, where Pzams 
is given by (Nelson fc Eggleton, 2001) 



„ 0.19JV/ii+0.47Mf;^ 
-TZAMS ^ — 



1+1. ISA!".-' 



(3) 



The fourth and fifth columns are the ages at which Roche 
lobe overflow begins (tRLOp) and the binary comes into con- 
tact (fcontact) in our calculation. The next three columns 
show some system parameters at tcontact, i.e. the mass of 
the primary Mi, the mass of the secondary M2 and the 
orbital period Pcontact- The remaining columns present the 
evolutionary results of the mergers: the lifetime on the main 
sequence (tMs), the central hydrogen mass fraction of the 
merger as constructed (Xucc) and after adjustment (Xncm), 
surface abundances for the elements H (Xhs), He (Xhos) and 
the ratio of C/N at the surface ((C/iV)s). 

Figure. [1] shows the location of the mergers on a colour- 
magnitude diagram (CMD) when the central hydrogen mass 
fraction reaches its maximum, at which we consider that the 
merger returns to thermal equilibrium and begins normal 
evolution. In a real case, the thermal equilibrium point of the 
constructed models are probably not just at the maximum 
of central hydrogen mass, but the divergence should be very 
small and have no influence on the results. '.' and 'x' in the 
figure represent the mergers from AS and AR, respectively. 
Open circles show the possible BSs produced from our mod- 
els in an old open cluster M67 and the dashed one shows 
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zero-age main sequence (ZAMS) for {Z,Y) = (0.02,0.28), 
where Z is metallicity and Y is helium content. In figure [1] 
we see that most of the mergers are located in the main se- 
quence while a few are to the left of ZAMS. Surface chemical 
composition and central hydrogen content are both respon- 
sible for placing objects to the left of ZAMS. Generally the 
mergers have a larger helium content than 0.28 because the 
matter from the secondary, including some helium-rich mat- 
ter produced via nuclear reaction, homogeneously mixes in 
the envelope of the primary. They are thus bluer than stars 
of the same mass and of the same central hydrogen mass 
fraction with Y = 0.28. Furthermore, the mergers may be 
on the left of ZAMS if central H contents close to 0.7, initial 
H mass fraction for Pop I. 

Sandquist, Bolte & Hernquist (I1997|) argued that a 
fainter BS should have a less massive progenitor, and there- 
fore has a lower helium content in the core. This means that 
BSs should spend an increasing amount of time near the 
ZAMS with decreasing mass, which will lead to a definite 
concentration toward the blue side of the main sequence in 
the region of low mass BSs. In Fig. [5] we present the dis- 
tribution of central hydrogen mass fraction at its maximum 
A'Hcm with the mass of the merger. No evidence in the figure 
shows that lower mass mergers should have higher hydrogen 
mass fractions in the center (say, more close to ZAMS). In 
our models, there are at least two factors affecting the cen- 
tral hydrogen mass fraction of a merger, i.e. the evolutionary 
degree of the parent stars when the system comes into con- 
tact and the development of the central convective region 
of the merger before normal evolution. The former is rele- 
vant to the mass of the progenitor and their contact ages. 
If we simply assume that all BSs in a cluster are formed 
at a same time (which is obviously not true), larger (more 
massive) progenitors have less hydrogen content in the cen- 
ter, but their mergers may develop larger convective regions 
in the centers, involving more of the H-rich matter from 
around the nuclear region of the primary in the center of 
the merger. As a consequence, the enhancement of central 
hydrogen content might be comparable to, or larger than, 
the less massive mergers. On the other hand, if a fainter BS 
is produced earlier than a brighter one - the formation time 
of a BS is strongly dependent on some other system's initial 
parameters as well as its progenitors' mass. The fainter BS 
has then evolved for a relatively longer time and may have 
less hydrogen in the center than the larger one. Therefore 
the blue concentration with mass decreasing predicted by 
Sandquist, Bolte & Hernquist (|1997p may be not true in a 
real case. In fact, we have not found the concentration obser- 
vationally, whereas there is a certain width of BSs sequence. 

We obtained some models with very little central hydro- 
gen content, i.e Xhchi very close to zero. Seen in Fig. [51 all 
these objects are from the AS channel. Long-time RLOFs of 
their progenitors lead the primaries be near the termination 
of main sequence when the systems being contact. Detailed 
calculations for these mergers show that their main-sequence 
timescales are in order of 10^ yrs. Such unusually long life- 
times are relevant to the development of central convective 
core. Comparison to the primaries of the progenitors, 
the products have higher masses and will develop 
larger convective core masses in the following evo- 
lutions, leading central hydrogen increase. On the 
other hand, nuclear reactions will consume some hy- 
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Figure 1. The locations of the mergers as they reach their max- 
imum hydrogen mass fraction on a colour-magnitude diagram. 
' . ' and ' X ' represent the mergers from AS and AR, respectively. 
Open circles show the merge models from Table 1 and the trian- 
gles are the objects with little central hydrogen content. 



drogen in the core. The maximum of central hydro- 
gen mass fraction is the equibibrium point at which 
the hydrogen involved in the core is comparable to 
that exhausted by nuclear reactions. After that, the 
consumed hydrogen is more than that involved in 
the core, but the growth of convective core adds 
fresh fuels and extends the lifetime in this phase. 
For a 1.6Mq star with a stable convective core, it is just 
about 2 X lO'^ yr when hydrogen mass fraction ranges from 
0.005 to 0.0000, while about 10* yr in these models. 

We should pay more attention to the mergers with little 
central hydrogen content. The triangles in figure 1 show the 
positions of these objects. From Table Al, we see that the 
primaries are very close to the exhausted of hydrogen. So it 
is very likely that the primaries have left the main sequence 
as eventually coalescence, and the mergers are subgiants and 
much redder than that shown in Fig. 1. Meanwhile, the de- 
velopment of convective core and the merger process are 
probably synchronous. After coalescence, the convective core 
develops more quickly than that in our models. Sandquist, 
Bolte & Hernquist (1997) once suggested that low-mass BSs 
(near the turnoff of a cluster) in M67 might be subgiant 
stars and a very small central hydrogen content in these 
objects may remarkably shorten their MS lifetimes and in- 
crease the chance of them being in the subgiant region. This 
may explain the observed spread in colour of low-mass BSs. 
These mergers with little central hydrogen in our models 
are similar to the objects mentioned by Sandquist, Bolte & 
Hernquist (|1997p . but they may be up to 2.7Mq, not just 
staying in the region around the turnoff in M67. 

Figures[3]and|3]show how the magnitude in the V band. 
My, and the colour, B — V , change with the merger mass 
Mbs at XHcm. My is calculated by 



Mv = Mboi - BC, Mboi = 4.75 - 2.5 x log(L/L0), 



(4) 



where BC and B — V are obtained 
by linear interpolation from the BaSel- 
2.0 model (|Lejeune, Cuisinier fc Buesel 1997] 



[Lejeune, Cuisinier fc Buesel 1998[ ). 
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Figure 2. The maximum of central hydrogen mass fraction for 
all the merger models. 



Figure 3. Visual magnitude changes with the mass of the merger. 
The dots and open crosses are for AS and AR, respectively. The 
stars represent the ones which are ruled out in eqs.(7) and (8). 



As seen in the figures, Mv and B — V are affected not 
only by the total mass, but also by the initial orbital period. 
For logPi/PzAMS =0.1, Figs. [5] and |6] present Mv vs Mbs 
and B — V YS Mbs, respectively. The solid lines are fitted 
by (two models marked with star are ruled out because they 
are outliers): 



-10. 93+22. 33M-2.65JW 



(P - V)o 



-10+10. 76Af 

_ 47.89+9.75Af-46.99M°'' 
-10+16. 07M 



(5) 
(6) 



The maximum errors for the two equations are 0.19 and 
0.019, respectively. 

Considering the initial orbital period, we get 



. 15. 87+19. 58A/-3.01A/^-' 



- 10+16. 22A/ 



-1.5(logPi/PzAMs-0.05) (7) 



(P - V)o = 



102+21. 73 Af- 101. 5 Af" 
-10+27. 44 



. 623(logPi/PzAMS-0-05) 
~^ 10+62^4a73 



(8) 



The maximum errors for the two equations are 0.29 and 
0.037, respectively. The four models marked with stars in 
FigsO and |4] are also ruled out in the fitting for the same 
reason as above. The distribution of error for eqs.(6) and 
(7), from detailed calculation, are shown in Figs[7]and[8] 

Initial parameters of the four systems ruled out in 
eqs.(6) and (7) are (Afi;, logg;, log(Po/PzAMs))= (0.89, 0.2, 
0.1), (0.89, 0.25, 0.1), (1.26, 0.15, 0.15) and (1.41, 0.25, 0..2). 
From Table Al, we see that all four of these mergers have 
central hydrogen content close to zero, i.e. 0.0007, 0.0020, 
0.0028 and 0.0007, respectively. This might be the reason 
that they deviate from most of our models 

The evolutionary tracks of the mergers with A/ — 
1.97M0 ((Mii,loggi)=(1.41, 0.40)) are shown in Fig.[9l Two 
stars with Y = 0.28 and M = 2.00 and 1.97A/0 are also pre- 
sented in the figure (dotted lines) for convenient comparison. 
In the figure, we can hardly distinguish the difference from 
various initial orbital period except for the main-sequence 
life, tMS, of the mergers. Though tus cannot be directly ob- 
served, it may affect observations, i.e. together with t 

contact 

of the progenitor system, it is crucial for whether or not the 
merger can be observed or not by now. At the same time, 
various initial orbital periods result in different luminosity 
functions for mergers - a longer initial orbital period leads a 




M/M^ 

Figure 4. Colour index B — V changes with the mass of the 
merger. The dots and open crosses are for AS and AR, respec- 
tively. 



shorter tms, but during the shorter tms, the merger mainly 
stays in the relatively high luminosity region (Fig. |9} . 



4 THE MERGERS AND BSS IN OPEN 
CLUSTERS 

If a BS in a cluster is the merger remnant of binary coa- 
lescence Q, then there are some constrains on the progeni- 
tor and on the merger itself. Firstly the total mass of the 
progenitor should be larger than the turnoff of the clus- 
ter Mto- Secondly, the contact age of the progenitor icontact 
should be appropriate. It cannot be so short that the merger 
has terminated its evolution on the main sequence, indicat- 
ing tcontact + tec + tus > ^cluster • Here ice IS the timescale 
from binary contact to final merger, fcontact also cannot be 
too long to complete its merger process at tcontact, indicat- 
ing tcontact + tec < tciustcr- If & binary has already come 



It just means AS and AR evolutionary channels in this section 
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Figure 5. Visual magnitude changes with the mass of the merger 
as log(Pi/PzAMs) = 0.10 . The dots and open crosses are for AS 
and AR, respectively. The stars represent the ones which are ruled 
out in fitting eqs.{5). The solid line shows the fitting curve via 
equation(5). 




Figure 6. Colour index B ~ V changes with the mass of the 
merger as log(Pi/PzAMs) = 0.10. The dots and open crosses are 
for AS and AR, respectively. The stars represent the ones which 
are ruled out in fitting eqs.{6). The solid line shows the fitting 
curve via equation(6). 



into contact but has not completed the mergering process 
(i.e. a W UMa system), it can also be considered as a BS 
dStryker 19931 ), e.g. S1036 and S1282 in M67. The initial 
parameter space for W UMa systems may be obtained by 

^contact ~t~ ^cc ^ ^cluster ^ ^contact . 

Figure [To] presents the lives of the mergers on the main 
sequence. We see that some low-mass BSs (i. e. M < 2.0Mq) 
may exist for about 10^ yr, which is long enough to be ob- 
served. In most cases, it is in order of 10* yr, which is similar 
to that of W UMa stars from observations, and therefore we 
may roughly estimate the contribution to BSs from AS and 
AR via the number of W UMa systems in a cluster. The 
estimation, however, is not absolutely since both of the two 
timescales have wide ranges and large uncertainties, and we 
cannot rule out other methods for creating W UMa 
systems except for AS and AR. Table 1 presents char- 
acteristics of some open clusters and the numbers of BSs 
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Figure 7. Visual magnitude error from equation (6) and the 
calculations. 
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Figure 8. Colour index B — V error from equation (7) and the 
calculations. 




Figure 9. Evolutionary tracks of the mergers after adjustment. 
(Afii, loggi)=(1.41, 0.40) for the solid lines and a, b, c, d, e repre- 
sent logPi/PzAMS = 0.10,0.15,0.20,0.25,0.30, respectively. The 
dotted lines show the evolution of two stars with surface helium 
content y=0.28 for = 1.97 and 2.OOM0. 
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are concentrated on low mass. Since the individual compo- 
nents almost have not evolved before contact, their mergers 
from this way have much longer timescales on the main se- 
quences. 



o 



LO 




o - 



1 1.5 2 2.5 3 3.5 

Figure 10. The timescale of the mergers staying on the main 
sequences. As tec = 5 X 10* yr, the open circles show the models 
fulfilled tcontact+ ice < 4.5x10^ yr and tcontact-l-tcc-l-tMS > 3.5x 
10^ yr. The squares and stars are for tcontact +tcc < 2.4 x 10^ yr 
and tcontact-l-icc+tMS > 1.4x10'-' yr and tcontact-l-tcc < 1.0x10^ 
yr and tcontact + tec + <MS > 0.9 X 10^ yr, respectively. 



Table 1. Characteristics of some open clusters (i.e. the age t and 
metallicity Z) and the numbers of BSs (A^bs), W UMa systems 
(A^wuMa) 3'Ud stars on the main sequence to two magnitudes be- 
low the turnoff {N2) in them. Most of the information is from 
Rucinski(1998) for A^UMai Ahumada & Lapasset (2007) for Abs 
and A2, and Xin & Deng{2005) for t and Z. The data not from 
the references above are marked with 1, 2, where 1 - Mochejska et 
al. (2004); 2 - Kafka et al.(2004). The stars means that we have 
not found related reports. 



ID 


t{Gyr) 


Z 




AwUMa 


N2 


Be 33 


0.7 


0.005 


2 


1 


270 


Tom 2 


1 


0.009 


17 


4 


440 


NGC 2243 


1.1 


0.007 


9 


2 


120 


NGC 2158 


1.2 


0.006 


40 


61 


700 


NGC 2660 


1.2 


0.02 


8 


0* 


150 


NGC 6939 


1.6 


0.02 


5 


62 


180 


NGC 3680 


1.6 


0.026 


1 


0* 


30 


NGC 752 


1.7 


0.014 


1 


1 


25 


NGC 7789 


2.0 


0.016 


22 


5 


130 


NGC 2682 


4.0 


0.02 


30 


3 


175 


NGC 188 


7.0 


0.024 


24 


7 


185 


NGC 6791 


7.2 


0.039 


75 


4 


800 


Be 39 


8.0 


0.01 


43 


9 


600 



and W UMa systems in them from observations. We see 
that A^'bs >> A'^wuMa in old clusters (i.e t > 2.0Gyr) and 
A'bs ~ A^wuMa in half of the left clusters, indicating that 
our models (binary coalescence from AS and AR) are not 
important for the produce of BSs in old open clusters, while 
likely play a critical role in some younger open clusters. In 
old open clusters, where stellar collisions may be ignored 
because of low stellar density, AML of low mass binaries 
is possibly dominated in producing BSs, since the time is 
long enough for binaries with initial orbital period about 2 
d evolving from detached to contact by AML and the merg- 
ers may be more massive than the turnoff. Meanwhile, from 
initial mass functions which have been presented, most stars 



4.1 Binary Samples 

To investigate BSs resulting from binary coalescence, we 
have performed a Monte Carlo simulation where a sample 
of 10^ binaries are considered (very wide binaries are actu- 
ally single stars) including BSs originated from AS and AR 
evolution channels. A single starburst is assumed in the sim- 
ulation, i.e. all the stars have the same age and metallicity 
{Z = 0.02). The initial mass function (IMF) of the primary, 
the initial mass ratio distribution and the distribution of 
initial orbital separation are as follows: 

i) the IMF of Miller & Scalo (IT979)) is used and the 
primary mass is generated from the formula of Eggleton, 
Fitchett & Tout (IT989|l : 

»,r 0.19X /q\ 

1 (l-X)0-''5-|-0.032{l-X)l/4 y'^! 

where X is a random number uniformly distributed between 
and 1. The mass ranges from 0.1 to lOOM©. 

ii) the mass ratio distribution is quite controversial and, 
for simplify, we only consider a constant mass ratio distri- 
bution (|Mazeh et al.l992|l . 

= 1,0 <(?< 1 (10) 

where q = M2/M1. 

iii) We assume that all stars are members of binary sys- 
tems and the distribution of separations is constant in loga 
(a is separation). 

/ N f ascp(a/ao)'", a < ao , , 

an(a) — < (11) 
[ Osop, ao < a < ai 

where a = 0.070, ao = 10i?Q,ai = 5.75 x W'^Rq = 0.13pc 
and m — 1.2. This distribution gives an equal number of 
wide binary systems per logarithmic interval and 50 per cent 
of systems with orbital periods of less than 100 yr. 

4.2 NGC 2682 

Some studies show that the metallicity of NGC 2682 
is a little different from the solar ([Carraro et al. 19961 
IFan et al. l"996|) . while some other studies con- 
cur that it is virtually indistinguishable from solar 
(|Hobbs fc Thorburn 19911 [FViel & Janes 1993J. So we select 
this cluster as the first sample to examine our models. 
The distance modulus m — M — 9.55 and the reddening 
E[B ~V) = 0.022 (|Carraro et al. 19961 [Fan et al. 199"6)l as 
we translate the theory results to observations. 

Previous studies on the age of this cluster showed 
several discrepancies. It may range from 3.2+0.4 Gyr 
(|Bonatto fc Bica 2003|l to 6.0 Gyr ( [Janes fc Phelps 1994| ). 
The study of VandenBerg fc Stetson (|2004|) derived an 
age of 4.0 Gyr. In the N-body model of this cluster 
( [Hurley et al. 2005] ), the authors investigated the behaviour 
around 4 Gyr. As mentioned in section 2, it is also a great 
uncertainty for the timescale tec . Table 2 and Table 3 present 
initial parameter spaces from difi'erent considerations for 
both BSs and W UMa systems in our grid. In Table 2, we 
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Table 2. Initial parameters for the mergers and W 
UMa systems in the grid for the old open cluster M67. 
We take tclustcr = 3.84 X lO^yr, corresponding turn-off 
is I.26M0. 





t^^ = 1.0 X 10* 


yr 








JVIii 






logPi 


the mergers 


1.26 




0.15 


0.1 


W UMa 


1.26 




0.05 


0.2 




Ccc — 0.1) X iU 


yr 






the mergers 


1.26 




0.25 


0.05 




1.41 




0.15 


0.2 


W UMa 


1.26 




0.05 


0.2 




1.26 




0.15 


0.1 




= 1.0 X 10® 


y^ 






the mergers 


1.12 




0.05 


0.05 




1.12 




0.3 


0.05 




1.26 




0.05 


0.15 




1.26 




0.1 


0.1 




1.26 




0.25 


0.05 




1.26 




0.35 


0.15 




1.26 




0.4 


0.15 




1.41 




0.2 


0.15 




1.41 




0.3 


0.15 


W UMa 


1.26 




0.05 


0.2 




1.26 




0.15 


0.1 




1.41 




0.15 


0.2 




1.41 




0.25 


0.15 



fixed the cluster age at tciustor = 3.8 x 10^ yr (corresponding 
to tiie turnoff of I.26M0) and varied tec, i-e tec ~ i x 10* 
yr, tec = 5 X 10* yr and tec = 1 x 10® yr. In Table 3, we set 
tec = 5 X 10* yr while the age of the cluster has a width from 
3.2 to 4.2 Gyr. then the conditions for valid parameter space 
to form BSs via case AS or AR are tcontact + tec < 4.2 x 10® 
yr, tcontact + tec + tus > 3.2 x 10® yr, and the constrains for 
W UMa systems in the cluster are tcontact < 4.2 x 10® yr 
and tcontact + tec > 3.2 X 10® yr. 

In Table [21 we see that tec will remarkably affect on 
initial parameter space for both of the two kinds of objects. 
For example, the initial parameter space is much larger for 
tec = 1 X 10® yr than those with other tec- We should bear in 
mind that, however, the long tec may be unreasonable from 
both of observations and theories (see section 2). The long 
tec may delay the appearance of the mergers and shorten 
their timescales on the main sequence, since individual evo- 
lution during mergering process are ignored. Especially for 
the primaries with a very littler hydrogen in the center at 
the system contact, the mergers may have never been on the 
main sequence. 

From TableO we see that, because of the existence of an 
age range, some models might be either mergers or W UMa 
systems. For large uncertainties of the age of the cluster, 
the initial parameter defined this way include almost all the 
models in Table [2] 

Some models with Mu — 1.41M0 appear in Table (2] 



Table 3. Initial parameters for the mergers and W 
UMa systems in the grid for the old open cluster M67 
(t = 3.2 - 4.2 Gyr and tee = 5 X 10* yr). 





the mergers 






W UMa 






log<3i 


logPi 




log9i 


logPi 


1.12 


0.05 


0.05 


1.12 


0.1 


0.05 


1.12 


0.3 


0.05 








1.26 


0.05 


0.15 


1.26 


0.05 


0.15 


1.26 


0.1 


0.1 


1.26 


0.05 


0.2 


1.26 


0.15 


0.1 


1.26 


0.1 


0.1 


1.26 


0.25 


0.05 


1.26 


0.15 


0.1 


1.26 


0.35 


0.15 


1.26 


0.3 


0.1 


1.26 


0.4 


0.15 








1.41 


0.05 


0.25 


1.41 


0.15 


0.2 


1.41 


0.15 


0.2 


1.41 


0.2 


0.15 


1.41 


0.2 


0.15 


1.41 


0.25 


0.15 


1.41 


0.25 


0.15 


1.41 


0.25 


0.2 


1.41 


0.3 


0.15 


1.41 


0.3 


0.2 


1.41 


0.35 


0.3 








1.41 


0.4 


0.3 









We checked the evolutionary details of them in Table Al 
and find that the system (A/n, loggi, log(Po/f'zAMs))= (1.41, 
0.15, 0.2) begins RLOF at 1.4 Gyr and becomes contact at 
3.3 Gyr. As tec = 5x 10* yrs, the merger is formed at 3.8 Gyr 
with Xucm = 0.0003 after self- adjustment and will leave the 
main sequence after 0.45 Gyr. As discussed in section 3, the 
merger of this system is very likely to be a subgiant when 
it formed. For other binaries with Mu — I.4IM0 in Table 
[21 though the primary is not near to zero, relatively long 
contact timescales (1 x 10®) lead the mergers be formed at 
about 3.5 to 3.8 Gyr. So they are still on the main sequence 
at the cluster age. 

Now we consider another case which is related to the 
BSs in a binary BS-I-MS in this cluster. It means that the sec- 
ondary is a BS before the system being contact. The parent 
binary for this case should fulfill the conditions as follows; 

(a) tcontact > tclustor > iRLOP; (b) at tdustor, M2 > Mto- 

Only one model is found in our results to possibly produce 
BS-I-MS in M67. The initial parameters for the model are 
(A//ii,loggi,logPi) = (1.41,0.25,0.2). At the age 3.82 x 10® 
yr, the secondary is 1.26A'/0 and the system comes into con- 
tact at 4x 10® yr with an orbital period of 0.5276 d. It means 
that the BS-I-MS may exist at least for 2 x 10*yr. A system 
like this should be near the turn-off if only the secondary 
(BS) is considered. However no observed BSs in this region 
are located in a binary like this. The evolved component 
may contribute to some luminosity l|Tian et al. 2006]) . 

By interpolating from Table Al, we get different num- 
bers of BS (A^Bs) under various conditions (Tabled). The 
location of the produced BSs on a CMD are presented in 

As mentioned before, AML of low mass binaries proba- 
bly makes a major contribution to BSs in this old open clus- 
ter. For low mass binaries, the individual components have 
almost not evolved (very close to ZAMS) before contact, and 
therefore their mergers may be replaced with ZAMS mod- 
els with a mass of the parent binaries. After coalescence, 
however, the mergers are much more massive and their evo- 
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Figure 11. Monte-Carlo simulation results for M67 at different conditions (Table |4]) 'x' represents the models from grid calculation, 
open circles show the models in Tableland the dots give the possible BSs. The observed BSs from Sandquist & Shetrone (2003) are 
also plotted in the figure. 



lutions cannot be negected agian. To examine the effect of 
AML in this cluster, we simply deal with low-mass binaries 
(Mii < IMq) as follows: 

(1) From the binary sample, we found out some binaries 
with Afii > l/2Mto and with the orbital period less than 0.5 
d at the cluster age by AML. The lower mass limit l/2Afto is 
sure that the mergers are probably larger than the turnoff, 
and the constrain on the orbital period leads the binaries 
possibly be contact before or at the cluster age. A semi- 
empirical formula for the orbital period variation is adopted 
here ( |Stepien 20061 ): 

fi^^ _(2.6 + 1.3) X 10-^0 p-i/3^-o.2P^,^ (^2) 
where Porb in days and time in years. For very short orbital 



periods the exponential factor is close to unity and varies 
very little during the subsequent evolution of the orbital 
period of the binaries. So it is ignored in this paper. 

(2) Find out the time, tRLOF, at which the primary is 
full of its Roche lobe for the selected binaries from step 1. A 
rapid binary evolution code written by Hurley, Pols & Tout 
(2002) is employed here and the Roche lobe of the primary 
J?cri is calculated by ( [Eggleton 1983[ ) 

Rcrl/A =0 692/°3+^n(l + ,l/3) ' (13) 

where A is the separation and q — M1/M2 

(3) Evolve the mergers to the cluster age. The mergers 
are ZAMS models instead with a mass of the systems. The 
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Table 4. Consequences of Monte Carlo simulation at 
different conditions for binaries undergoing AS and 
AR evolutionary channels in NGC 2682 (M67). 







conditions 




A^BS 




case 


1 


tec = 1 X 10*yr 
tclustor = 4 X 10^ yr 




5 


3 


case 


2 


tec = 5 X 10*yr 
tciustor = 4 X 10^ yr 




5 


10 


case 


3 


ice = 1 X lO^yr 
tciustor = 4 X 10^ yr 




14 


15 


case 


4 


ice = 1 X 10*yr 
tclustor = 3.4-4.2 X 10^ 


yr 


20 


13 


case 


5 


tc^ = ^x 10*yr 
tclustor = 3.4-4.2 X 10^ 


yr 


25 


19 


case 


6 


ice = 1 X lO^yr 
tclustor = 3.4-4.2 X 10^ 


yr 


38 


36 



Table 5. AML Consequences for NGC 2682. 







conditions 


A^BS 


WwUMa 


case 


1 


tec = 1 X 10**yr 
iciuater = 4 X lO" yr 


855 


14 


case 


2 


tec = 5 X lO^yr 
iclustcr = 4 X 109 yr 


964 


84 


case 


3 


ice = 1 X lO'-'yr 
tclustcr = 4 X 109 yr 


1079 


204 


case 


4 


ice = X lO^yr 
iclustor = 4 X 109 yr 


824 





case 


5 


tj.^ = 5 X lO^yr 
^cluster = 3.5 X 10^ yr 


996 


120 


case 


6 


ice = 5 X lO^yr 
tclustor = 4.5 X 10^ yr 


1194 


109 



starting age is iaLOF, since the systems will reach contact 
very quickly from the onset of RLOF because of AML. Dif- 
ferent timescales tec are adopted here. Some of the binaries 
from step 1 have left the main sequence, but most of them 
remain on the main sequence. 

Table 5 gives the results from different cluster age as- 
sumptions and different tec, and the magnitudes and colors 
for the mergers from these cases are also presented in Fig. 
11. 

In Fig llll we see that, the products from binary evolu- 
tion (AS and AR in the paper) and those from AML cover 
different regions. The mergers from conservative evolution 
are located in the region of high luminosity, i.e far away 
from the the turnoff, and the remnants from AML occupy 
the region with low luminosity and have a scatter on the 
color. The figure also indicates that AML of low-mass bi- 
naries is much more important in this old open cluster. In 
some cases (case 3 to 6), the produced BSs from conserva- 
tive evolutionary appears a certain width with decreasing 



mass. Another phenomena worth noticing is in case 6, the 
produced BSs result extends very close to F81, the observed 
brightest BSs in this cluster. Further study shows that the 
lower limit of the cluster age may be 3.8 Gyr for ice = 1 x 10^ 
yr to result in the produced BSs extending to the position 
of F81. 

As shown in Table 4, we only obtained 38 BSs from a 10® 
binary sample in the widest condition. NGC 2682 is an open 
cluster and has a much less stars than 10®. In Hurley et al. 
([2005)) . 12000 single stars and 12000 binaries (iV = 24000) 
are adopted in the best fitting model. Multipying this factor, 
we just got 1 BS with high luminosity. There are 3 BSs from 
the observations, however, located in the region covered by 
the mergers from AS and AR. The number of stars initially 
in this cluster is important here, but it is difficult for us to 
estimate it because of the incompleteness of observations as 
determing the current mass M and the uncertainties of mass 
loss history of the cluster when converting M to an initial 
mass Afo. IMF may also affect the results. We adopt the 
field single star IMF for the primaries and the lower mass 
limit is O.IM0 in the binary sample, however, M67 is an 
old open cluster and rich in binaries, and then it possibly 
has different IMF as well as different lower mass limit as 
described by Chabrier ([2003]! . 

In the region with low luminosity, i.e the region covered 
by the mergers of AML, the numbers shown in Table 5 are 
enough to account for BSs in this region. Dynamically un- 
stable mass transfer (AD) is also possible to form BSs in 
this region, since the binaries undergoing AD have larger 
initial mass ratios in general, leading the smaller mergers. 
According to the result of Nelson & Eggleton (2001), we 
roughly estimate the number of binaries undergoing AD evo- 
lution with initial primary mass between 0.89 and I.26M0 
as N = 10®. It is about 41-63 and the mass of the bina- 
ries are less than 1.36Afto, indicating that AD may produce 
BSs near the turnoff indeed, but the contribution is much 
smaller than that of AML in low-mass binaries. 

In Hurley et al. (l2005p . the authors obtained seven BSs 
from case A mass transfer leading to coalescence by un- 
perturbed evolution. The age of the mergers in that paper 
were calculated based on the assumption that core hydrogen 
burning proceeds uniformly and that the end of the main se- 
quence is reached when 10 per cent of the total hydrogen has 
been burnt. These assumptions for calculating the MS lives 
of the mergers are questionable as they mentioned in the pa- 
per. According to our calculations, almost all of the mergers 
have relatively short main-sequence lives, i. e. less than 10^ 
yr (see Fig llOp . It means that five of the seven BSs are likely 
to leave the main sequence and no longer be observed as BSs 
at 4 Gyr. The range of age for M67 is likely to improve the 
result. 

4.3 NGC 2660 

NGC 2660 is an intermediate-age open cluster with a solar 
metallicity. The early study by Hartwick & Hesser (1973) 
showed the following properties for this cluster: E{B — V) = 
0.38, (m - A/)o = 12.3 ± 0.3, age ~ 1.2 Gyr, metallicity 
similar to the Hyades, and high possibility of membership 
for the N-type carbon star. There are some uncertainties for 
the determination of the cluster age, e.g. Lynga (1987) cited 
1.6 Gyr, Jane & Phelps (1994) gave an age of 0.9 Gyr, while 
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Table 6. Consequences of Monte Carlo simulations from different assumptions for binaries undergoing AS and AR 
evolutionary channel in NGC 2660. 







conditions 




A^BS 

(Ml = O.8M0) 


A'wUMa 

(Ml = O.8M0) 


A^BS 

(Ml =O.1M0) 


A'wUMa 

(Ml =O.1M0) 


case 


1 


tclustcr = 1-2 X 10^, tec = 5 X 10'^ yr 




268 


59 


70 


11 


case 


2 


tclustcr = 1-2 X 10^ tec = 1 X 10* yr 




284 


117 


71 


21 


case 


3 


tclustor = 1-2 X 10^, ice = 5 X 10» yr 




296 


631 


72 


126 


case 


4 


<cluster = 1-2 X 10^ ice = 1 X lO^ yr 




87 


1038 


21 


220 


case 


5 


tcluster = 0.9 - 1.2 X 10^ ice = 5 X lO'^ 


yr 


693 


454 


141 


87 


case 


6 


tcluster = 0.9 - 1.2 X 10^, tec = 1 X 10* 


yr 


703 


525 


149 


106 


case 


7 


icluster = 0.9 - 1.2 X 10^ tec = 5 X 10* 


yr 


445 


934 


102 


194 


case 


8 


tclustcr = 0.9 - 1.2 X 10^, tec = 1 X 10^ 


yr 


87 


1125 


21 


241 



Carrio & Chiosi (1994) derived 0.7 Gyr. The latest report 
on this cluster is Sandrelli et al. (1999). It was shown that, 
metallicity about solar, (m-M)o = 12.1-12.3, E{B-V) = 
0.37 — 0.42, age < IGyr, with a fraction of binaries of about 
30 per cent. 

According to the new catalogue of blue stragglers in 
open clusters ( |Ahumada &: Lapasset 2007| , there are 8 BSs 
in this cluster while the cluster age is logt — 9.03 (about 
1.1 Gyr). The turnoff is around 2M0 for this age, indicating 
that AML in low-mass binaries has little influence on the 
product of BSs in this cluster. Table 6 presents our simu- 
lation results from different assumptions of the cluster age 
and tec by interpolating in Table Al. 

Though we show the simulation results for tec = 5x 10* 
and tec ~ i X 10^ yr in Table |6l we should bear in mind 
that NGC 2660 is an intermediate-age cluster and a long 
tec is unreasonable for it from both observations and our 
models. Such a long tec as 5 x 10* or 1 x 10^ yr will inevitably 
cause large differences of both components of a binary in our 
models, and therefore a shorter tee, namely, ice = 1 x 10* 
and tec = 5 x lO'^ yr, is more appropriate for this cluster. 
From Table [6] we see that BS number is not sensitive to tee 
in this cluster, but the contact systems are very sensitive to 
tec, since a longer tec indicates a larger parameter space for 
these systems. Contact binaries exist in the cluster for all the 
cases in Table[Sl but it is less than the BSs number produced 
from AS and AR when tee = 1 x 10* and tec = 5 x lO'^ yr, 
especially in case 1, it is just of 1/6 to 1/5 of that of BSs. 

There are 150 stars on the main sequence to two magni- 
tudes below the turnoff in NGC 2660 and 175 in NGC 2682, 
as shown in Ahumada & Lapasset l|2UU7|l . From this com- 
parison, our models may produce the necessary BSs in NGC 
2660 when Ml = 0.8Mq. If Ml = O.IMq for the binary 
sample, however, the produced BSs number from AS and 
AR is just 1/5 to 1/4 of that from Ml = O.8M0. Normal- 
ization is necessary here to estimate the BSs birthrate from 
our models in this cluster and the same problem mentioned 
in section 4.2 appears here again. Nevertheless, our models 
may explain several BSs in this intermediate age cluster. 

Figure [T^] shows some examples from Monte 
Carlo simulation on CMD. The observed BSs 
( |Ahumada fc Lapasset 2007[ ) and an isochrone {t = 1.2Gyr) 



are also plotted in the figure. Several BSs are located in the 
region of our models from AS and AR. In the region with 
the most chance from our models Mv « 1.0, however, there 
are no BSs observed at present. The observed peak appears 
about Afv ~ 2.0, one magnitude lower than the theoretical 
value. Mass loss during merger process may shorten the 
discrepancies, but we should find out an appropriate 
mechanism to explain such a large mass loss about O.5M0, 
which is much larger than that from observations and from 
smooth particle hydrodynamic simulations. 

Though the BSs number (or birthrate) from our models 
is not sensitive to tec , the initial parameter space will change 
with tee, resulting in some differences for the mergers. In the 
upper right panel, we show the results from case 1 to case 
3 with Ml = O.IMq. As seen in this figure, the mergers 
from a long tee, i.e. 5 x 10* yr, extend to a higher luminosity 
and concentrated on the ZAMS, but small differences appear 
between tec = 1 x 10* yr and tec = 5 x 10^ yr. 

We have not considered the effect of AML in 
this cluster, since we simply assume that the con- 
servation assumptions are reasonable for stars with 
spectra from GO (about l.OSM©) to Bl according 
to Nelson &: Eggleton (|200ip . In fact, a main se- 
quence star can generate a magnetic field as long 
as it has a sufficiently thick convective envelop, i.e. 
with a spectrum later than F8. In the study of An- 
dronov, Pinsonneault & Terndrup (|2006|) . the mass 
threshold for AML is between 1.2 and I.4M0, then 
the total mass of the merger might be up to 2.4 
and 2.8A/0, above the turnoff of this cluster. So it 
could be possible that AML is also impacting the 
blue straggler population in this intermediate-age 
cluster. This provides another possible explanation 
for the magnitude offset from the models and obser- 
vations. 



5 DISCUSSIONS AND CONCLUSIONS 

In Sect. 4, we notice that the timescale from contact to com- 
plete coalescence, tee, strongly affects the initial parameter 
space of primordial binaries which eventually produce sin- 
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Figure 12. Monte-Carlo simulation results for NGC 2660 from different assumptions (Table[6ll. The dots show the blue stragglers listed 
by Ahumada & Lapasset (2007), and the error bars are from Sandrelli et a. (1999), i.e. m - M = 12.1 - 12.3, E(B -V) = 0.37 - 0.42. 
The stars present the positions of BSs when m — M = 12.2, E{B — V) = 0.31 as shown in Ahumada & Lapasset (2007). In the upper 
right panel, the circles, dots and triangles are for tec equal to 5 X 10^, 1 X 10* and 5 X 10®, respectively. 



gie BSs in a cluster. On the other hand, there are some 
conflicting estimates for tec from observations and theo- 
retical models. In this paper we adopt empirical val- 
ues, i.e. tec is short in comparison to the evolution 
timescale of both components in a binary, and ig- 
nore the changes during the merger process. In this 
section, we will first discuss the consequences of a 
long tec- Because of evolution of both components 
during the merger process, the primaries have lower 
central hydrogen content and the matter from the 
secondaries have larger He content. The former re- 
sults in a redder colour for the mergers while the 
latter makes the mergers bluer. So the final posi- 
tions of the mergers are possibly similar to those 
shown in this paper, except that the primaries have 
left the main sequence at final coalescence. This case 
will appear in the mergers with little central hydro- 
gen. For example, a star with 2A/q may evolve from 
ZAMS to exhausted of central hydrogen in 10^ yr, 
and then none of the mergers from binaries with pri- 
mary' masses larger than 2Mq will be on the main 
sequence if tec = 1 x 10^ yr. For the primaries with 
very little hydrogen in the center at contact, the 
mergers may never be on the main sequence even 



in the cases of short tec- The long tec also delays 
the appearance of the mergers and shortens their 
timescales on the main sequence. The latter has not 
been exactly expressed in our models, and therefore 
we just see that the mergers from a long tec have 
larger luminosities as shown in section 4. 

In our binary evolutions, we have not included 
AML, which exists in low-mass binaries and may be 
the main course making the binaries change from de- 
tached to contact and finally coalesce, resulting in a 
large contribution to BSs in old clusters, e.g NGC 
2682. In young and intermediate-age clusters, how- 
ever, AML has little contribution to the birthrate of 
BSs, since (a) the time is not long enough for bina- 
ries to go from detached to contact and (b) the mass 
of the mergers is probably less than the turnoff of 
the cluster even though their parents may coalesce 
in the cluster age. So, we simply estimated the ef- 
fect of AML in NGC 2682 while negecting it in NGC 
2660. 

The mass loss during the merger process can also 
affect our result, mainly the location on the CMD of 
the products. As shown in NGC 2660, no BSs have 
been observed in the region with the most opportu- 
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nity from our models. Because of mass loss, the merg- 
ers will be fainter than those given in the paper. How- 
ever the faintness will be slight since the mass loss is 
not vast from both observations and smooth particle hy- 



drodynamic simulations ( 
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carry some angular momentum out from the parent bi- 
nary. By analyzing the BS spectra from Hubble Space Tele- 
scope (there is an apparent continuum deficit on the short- 
wavelength side of Balmer discontinuity ), De Marco et al. 
H2004P argued that some BSs might be surrounded by a cir- 
cumstellar disk. However, Porter & Townsend (2005) showed 
that the flux deficits may be attributed wholly to rapid ro- 
tation. The rotation rates needed are of the order of those 
found in the study of De Macro et al. H2005p . Whether the 
flux deficits shortward of the Balmer jump are induced by 
a circumstellar disk or rapid rotation, it provides a possible 
explanation for the orbital angular momentum of the system 
after coalescence. Such a large mass loss as shown in NGC 
2660 (about O.SM©), however, is a problem and should be 
explained reasonably in physics. 

Based on some assumptions, we studied the mergers of 
close binaries from AS and AR evolution by detailed evo- 
lutionary calculations. The products from our models may 
stay on the left of the ZAMS and have no central concentra- 
tion with decreasing mass. Because of the development of 
the convective core, the mergers with little central hydrogen 
(less than 0.01) in our models have unusually long timescales 
on the main sequence (~ 10* yrs). These objects are proba- 
bly subgiants as they are formed, since the primaries in the 
progenitors also have little central hydrogen and may have 
left the main sequence during merger process. 

The mergers from our models stay on the main sequence 
for a timescale in order of 10* yrs. Some low-mass mergers 
may stay on the MS for about 10^ yrs. The timescale is simi- 
lar to that of W UMa stars from observations, and therefore 
we may roughly estimate the contribution to BSs from AS 
and AR via the number of W UMa systems in a cluster. The 
estimation, however, is not absolutely since both of the two 
timescales have wide ranges and large uncertainties, and we 
cannot rule out other methods for creating W UMa systems 
except for AS and AR. Comparison to observations indi- 
cates that our models (binary coalescence from AS and AR) 
are not important for the produce of BSs in old open clus- 
ters, while likely play a critical role in some younger open 
clusters. 

We performed Monte Carlo simulations to examine 
our models in an old open cluster NGC 2682 and in an 
intermediate-age cluster NGC 2660. The effect of AML was 
estimated in NGC 2682 in a simple way, where the merg- 
ers are replaced with ZAMS models. In NGC 2682, binary 
mergers from our models cover the region with high lumi- 
nosity and those from AML are located in the region with 
low luminosity, existing a certain width. The BSs from AML 
are much more than those from our models, indicating that 
AML of low mass binaries makes a major contribution to 
BSs in this cluster. Our models are corresponding for sev- 
eral BSs in NGC 2660. In the region with the most oppor- 
tunity on CMD, however, no BSs have been observed. Our 
results are well-matched to the observations if there 
are ~ 0.5Mq of mass loss in the merger process, but 



a physical mechanism for this much mass loss is a 
problem. 
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sequence {tMs), the central hydrogen mass fraction of the 
merger as constructed (Xhcc) and after adjustment (Xncm), 
surface abundances for the elements H (Xhs), He (Xhos) and 
the ratio of C/N at the surface ((C/iV)s). 



APPENDIX A: EVOLUTIONARY RESULTS IN 
THE PAPER 

We choose a set of binaries undergoing AS and AR evolu- 
tion from Nelson & Eggleton (|200ip to study the charac- 
teristics of the merger products and their connections with 
blue stragglers. Table Al gives the initial parameters of 
the binary systems, their RLOF information and the struc- 
tures and evolutionary consequences of the mergers. The 
first three columns contain the initial mass of the primary 
Mii, the initial mass ratio (the primary to the secondary) 
and the initial orbital period Pi in logarithmic, where Pzams 
is given by (Nelson & Eggleton, 2001) 



ZAMS 



1+1. 18M-^ 



(Al) 



The fourth and fifth columns are the ages at which Roche 
lobe overflow begins (faLOp) and the binary comes into con- 
tact (tcontact) in our calculation. The next three columns 
show some system parameters at tcontact, i.e. the mass of 
the primary Mi, the mass of the secondary M2 and the 
orbital period Pcontact- The remaining columns present the 
evolutionary results of the mergers: the lifetime on the main 



